Quantum Quench in the Transverse Field Ising Chain 
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We consider the time evolution of observables in the transverse field Ising chain (TFIC) after a 
sudden quench of the magnetic field. We provide exact analytical results for the asymptotic time and 
distance dependence of one- and two-point correlation functions of the order parameter. We employ 
two complementary approaches based on asymptotic evaluations of determinants and form-factor 
sums. We prove that the stationary value of the two-point correlation function is not thermal, but 
can be described by a generalized Gibbs ensemble (GGE). The approach to the stationary state can 
also be understood in terms of a GGE. We present a conjecture on how these results generalize to 
particular quenches in other integrable models. 

PACS numbers: 



Recent experiments on trapped cold atomic gases p] [2] 
have raised intriguing fundamental questions regarding 
the non-equilibrium dynamics of correlated many-body 
quantum systems. These cold atom systems are suffi- 
ciently weakly coupled to their environments as to al- 
low the observation of essentially unitary nonequilibrium 
time evolution on long time scales. The quantum New- 
ton's cradle experiments of Kinoshita et al [5] in partic- 
ular have focussed the attention on the roles played by 
dimensionality and conservation laws. The observed ab- 
sence of "thermalization" in quasi one dimensional con- 
densates was attributed to the experimental system be- 
ing approximately describable by a quantum integrable 
many-body theory. This in turn initiated vigorous re- 
search on clarifying the role played by quantum integra- 
bility in determining the stationary (late time) behaviour 
of nonequilibrium evolution in correlated quantum sys- 
tems [4TllO| . The simplest way of driving a quantum sys- 
tem out of equilibrium is by means of a quantum quench: 
a system is prepared in the ground state of a given Hamil- 
tonian H(ho), where ho is an experimentally tuneable 
parameter such as a bulk magnetic field. At time t = 
the parameter ho is changed suddenly to a different value 
h and one then considers the unitary time evolution of 
the system by means of the new Hamiltonian H(h). Cen- 
tral issues that have been investigated are whether the 
system relaxes to a stationary state, and if it does, how 
to characterize its physical properties at late times. It 
is widely believed (see e.g. [3] for a comprehensive sum- 
mary) that the behaviour of local observables (such as 
one and two-point correlation functions) can be described 
in terms of either an effective thermal (Gibbs) distribu- 
tion or a GGE [4 . It has been argued that the latter 
arises for integrable models, while the former obtains for 
generic systems [1H7]. However, several recent studies [5]- 
ITU] suggest that the behaviour is more complicated and 
in particular depends on the initial state. Moreover, open 
questions remain even with regards to the very existence 
of stationary states. For example the order parameter of 



certain mean-field models have recently been shown to 
display persistent oscillations [11]. 

Two recent works have raised another crucial issue in 
the debate on thermalization, namely the role played by 
the considered observables [H] ■ More precisely, it was 
pointed out that the locality of the observable with re- 
spect to the elementary excitations is expected to affect 
the late time behaviour of an observable after a quan- 
tum quench. In light of the available experimental and 
theoretical results further clarification of the role of inte- 
grability on the time evolution after a quantum quench 
calls for exact analytical results on "generic" correlation 
functions, i.e. those corresponding to observables non- 
local with respect to the elementary excitations. 

In the following we present such results for the par- 
ticular case of the transverse field Ising chain, which 
is a crucial paradigm for quantum critical behaviour. 
While the model admits a representation in terms of 
free fermions, the order parameter is non-local with re- 
spect to the fermionic degrees of freedom, which ren- 
ders it an ideal testing ground for thermalization ideas. 
Although the model has been widely analyzed in the 
context of quantum q uenches 1 1 3H1 7] . the non-equilbrium 
evolution of order parameter correlation functions are 
still not known analytically. In this letter we present an- 
alytical results for the full asymptotic time and distance 
dependence of one- and two-point correlation functions 
of the order parameter in the thermodynamic limit af- 
ter a quantum quench within the ferromagnetic phase. 
We also present partial results for quenches within the 
paramagnetic phase and across the critical point. Our 
results are obtained by two independent, novel methods. 
The first is based on the determinant representation of 
correlation functions characteristic of free-fermionic theo- 
ries. The second is based on the form-factor approach [T5] 
and is applicable more generally to integrable quenches 
in interacting quantum field theories [19] . This method 
complements existing analytical/semi-numerical meth- 
ods used for studying quantum quenches in integrable 
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systems [JJ HU] , but has the advantage of providing ana- 
lytic answers directly in the thermodynamic limit. 

The model. We consider the spin-^ TFIC Hamiltonian 
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The Hamiltonian can be diagonalized by a combination 
of Jordan- Wigner and Bogoliubov transformations (see 
e.g. [15]). The dispersion of the elementary fermion 
excitations is e^k) — y/h 2 — 2h cos k + 1. The system 
is initially prepared in the ground state at a field h . 
The field is then instantaneously changed from ho to h 
and unitary time evolution with Hamiltonian H{h) en- 
sues. We are interested in the time evolution of the or- 
der parameter p x (t) = (af(t)) and its two-point function 
p xx (£,t) = (af (t)crf,g(t)) . Due to translational invari- 
ance the 1-point function is position independent and the 
2-point function depends only on the distance £. An im- 
portant role is played by the difference A& of the Bogoli- 
ubov angles diagonalizing H{h) and H(h ) respectively 
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Quenches within the ordered phase (h, h < 1). We 
find that at late times the order parameter relaxes to 
zero exponentially fast 
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where e' h (k) — deh(k)/dk. The two-point function of the 
order parameter exhibits exponential decay both in time 
and distance (8(x) denotes the Heaviside function) 
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In the £ — > oo limit the first factor is equal to unity and 
p xx (oo,t) = (p x (t)) , confirming cluster decomposition 
in our non-equilibrium situation. Fig. [T] shows a com- 
parison of our asymptotic result for p xx (£,t) to numeri- 
cal data, establishing the accuracy of the former even for 
relatively short separations and times. We note that Q 
holds even for quenches to or from the quantum critical 
point and agrees with the general form put forward in [3] 
on the basis of semiclassical arguments. 

The stationary state. The result Q allows us to make 
exact statements regarding thermalization in the model. 
The one-point function is trivially thermal, since it van- 
ishes for t — ¥ oo as was already pointed out in [5] I17j . 
On the other hand, in this limit the two-point function 




FIG. 1: p xx {£,t) for the quench ho = 0.3 -> h = 0.5 at fixed 
distance I — 20 and I = 40 against the prediction in Q. 
The overall amplitude of p xx has been used as the same fit 
parameter in both cases. 



exhibits exponential decay with a correlation length 
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This is reminiscent of the behaviour of the equilibrium 
two-point function p X q(£) at temperature T, which de- 
cays exponentially with correlation length (for h < 1)[T5] 
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For p xx (£, oo) to be thermal, £ would have to equal £r eff , 
where the effective temperature T e g is determined by the 
requirement that the (average) energy in the initial state 
(ipo\H(h)\tpo) is given by the thermal average (H(h))T e{{ - 
This leads to the following equation fixing T e g 
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With T off given by (M) we find that ^ and 

p xx (£, oo) is therefore never thermal. On the other hand, 
for small quenches corresponding to small values of T ef f , 
£ and £-r cff can be seen to coincide to order (h — h ) 2 
[24] . Hence there exists a small quench regime, where 
the thermal result provides a good approximation. This 
agrees with the numerical findings of Ref. [UJ, which 
were reported to be consistent with thermal behaviour 
at low effective temperatures. In [3] it was proposed that 
the stationary state of integrable models can be described 
in terms of a GGE, defined by maximizing the entropy 
while keeping the energy as well as higher conservation 
laws fixed. For the particular case of the TFIC this re- 
sults in a mode- dependent effective temperature given by 



cos At = tanh 



2T cff (fc) 
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Inserting this relation into the expression © for £ T in 
Eq. (|6| results in £ ([5]), which proves that the station- 
ary behaviour after a quench is indeed described by a 
GGE. To the best of our knowlege this demonstrates for 
the first time that the GGE applies even to correlation 
functions of observables that are nonlocal with respect to 
the elementary excitations in a lattice model, thus signif- 
icantly generalizing the results of [4]-[6|. Our expression 
([5]) for the correlation length in terms of T c g(k) ^ ex- 
plains nicely observations made in [17] , that numerical 
results on p xx are described more accurately by a mode- 
dependent correlation length than by the thermal result 
In [7 it was argued that the GGE quite generally 
describes the stationary behaviour of one point functions 
for quenches in integrable quantum field theories. Tak- 
ing the appropriate scaling limit of our result suggests 
the correctness of the assumptions made in [7] . 

Approach to the stationary state: Suprisingly the decay 
time in ^ can be similarly explained in terms of a GGE, 
even though it is not a property of the stationary state. 
The one point function ([3| is characterized by exponen- 
tial decay with rate = — J* ^e' h (k) In (cos A/.). This 
can be interpreted as the average mode-dependent decay 
time T _1 (fc) = e' h (k)( y ^ 1 (k), obtained by multiplying the 
mode-dependent inverse correlation length by the veloc- 
ity. The relaxational behaviour of the two-point function 
can be understood following [5] by rewriting Q as 



(p x (t)) 2 





' r dk 


r I 2t I 


exp 


Jo T 


k(k) r(k)\ 



9(2e' h (k)t-£) 



(9) 

The thcta-function expresses the fact that a given mode 
can only contribute to the relaxational behaviour if the 
distance £ lies within its forward "light cone", while the 
form of the remaining factor follows from the known sta- 
tionary behaviour. Numerical studies of the characteris- 
tic coherence time for the non equal-time two point func- 
tion were found to be compatible with thermal behaviour 
|17j . It would be interesting to revisit this analysis in 
light of our current findings. 

Quenches originating or ending in the disordered 
phase. Here the behaviour of correlation functions of 
o~ x is more involved |17j . A complete summary of our 
results on the dynamics in this case is beyond the scope 
of this note and will be reported elsewhere [19]. On the 
other hand, the correlation length characterizing the sta- 
tionary behaviour of p xx (£,t = 00) ~ exp(— £/£) for an 
arbitrary quench can be cast in the simple form 

C 1 = 0(h - l)9(h - 1) In (mm[h , hi}) 
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where x± — j[min(h, h 1 ) ± l][min(/io, h ) ± 1] and 
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This agrees with &4 1 and the 
and ho — 00 [TS]. Crucially, it 



can be proved by a direct calculation in the framework 
of the Toeplitz determinant approach summarized below 



that the result ( 10 1 agrees with the predictions of a GGE. 
We believe this calculation generalizes straightforwardly 
to nonlocal correlators in other free fermionic theories, 
which suggests that the GGE correctly predicts infinite 
time behaviour for both local and nonlocal observables 
in such theories. 

Method I: Determinant approach. We focus on the two- 
point function p xx (£, t) , which can be written as the de- 
terminant of a 2£x 21 block Toeplitz matrix T [l3 l fTfl fT9] . 
The matrix elements of T depend explicitly on the time 
t. In the stationary state the t dependence disappears 
[T5] and the large-£ behaviour can be obtained by ap- 
plication of the generalized Szego lemma |T3]. resulting 
in (10). The dynamics in the limit t,£ — > 00 at fixed 



ratio tjl is much more difficult to determine, as the el- 
ements of T then depend on the matrix dimension itself 
and Szego's lemma does not apply. To deal with this 
situation we employ a method similar to |16j . In order 
to calculate lndet|T| = Trln|T|, we consider the mo- 
ments of T, i.e. TrT 2 ™ (we find that odd moments are 
subleading). Calculating these moments gives 
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where s(x) — x6{x). The trace of any analytic function 
/ of T can be formally expanded in the moments. In our 
case we are interested in f[x) = In |x|, which is analytic 
in the principal strip for any x 7^ 0. As the symbol of the 
block Toeplitz matrix has winding number zero about 
the origin and cos(Afc) is always non zero we can resum 
the expansion of the logarithm to obtain Q. We have 
generalized these results to the case of the XY spin chain 
in a field [33]. 

Method II: Form-factor approach. This approach ap- 
plies more generally to quenches in integrable (interact- 
ing) quantum field theories. We focus on the 1-point 
function in the ordered phase. The ground state for 
\h\ < 1 spontaneously breaks the Z2 symmetry of the 
TFIC, resulting in an initial (ground) state of the form 
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where R and NS refer to the periodic/antiperiodic sec- 
tors of the free-fermion theory respectively and e.g. 
\B) K = exp(iJ2 0<peR K(p)blbl p ) |0) R , where K(p) = 

tan [Ap/2] and b' p is a fermion creation operator with mo- 
mentum p. The 1-point function is 
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Expanding the "boundary states" |-B)r,ns results in a 
Lehmann representation for (13 1. Crucially, the matrix 
elements (form factors) of cr^(i) between multifermion 
Hamiltonian eigenstates are known exactly for the TFIC 
[22j . The main idea for evaluating the Lehmann repre- 
sentation then follows the finite temperature case [2"T1|23| . 
For a small quench the (total) density uq of fermion exci- 
tations in the initial state constitutes a small parameter. 
In this case one can use the K(k)-m&trix as an expansion 
parameter. One then observes that the form factors ap- 
pearing in the Lehmann representation are singular when 
momenta in the in and out states coincide. The leading 
(in the density no) contribution to the one point fuction 
is obtained by summing all terms with the strongest sin- 
gularities at a given order in the expansion in powers 
of K. This amounts to the exponentiation of infrared 
singularities. We note that just as in the finite temper- 
ature case [23] infinite volume divergences encountered 



supported by the EPSRC under grant EP/D050952/1 
(FHLE) by the ESF network INSTANS (PC and MF). 



in evaluating the numerator of ( 13 ) cancel against anolo- 



gous divergencies in the denominator. The result of these 
calculations is 



oc exp 
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., i —K 2 (k)2e'(k) 



(14) 



The decay rate agrees with the the leading term in the 
expansion of ^ in powers of K 2 (k). The correction to 
the K 2 (k) factor in |l4| ) are found to be 0(K 6 ), again 
in agreement with ([3J. We note that ( [M] ) provides an 
excellent approximation to ([3| as long as h, ho are not 
too close to the critical point. The two-point function can 
be analyzed in an analogous manner |19j and the results 
again agree with the appropriate expansion of Q. 

Summary and Discussion. We have obtained exact an- 
alytic results for the long distance and time asymptotic 
behaviour of one and two point functions of the order 
parameter a x in the TFIC atfer a quantum quench. We 
have shown that the stationary expectation value of the 
two-point function is not thermal, but can be described 
by GGE. We have further shown that the approach to the 
stationary state for quenches within the ordered phase 
can be understood in terms of a GGE as well. Our work 
further demonstrates the importance of having analytic 
results at one's disposal when trying to draw conclusions 
regarding the statistical description of stationary state 
properties. Finally, we comment our newly developed 
method based on form factors generalizes directly to in- 
tegrable quenches in integrable quantum field theories. 
These are characterized by the requirement that the ini- 
tial state is compatible with factorizable scattering. We 
conjecture that the stationary behaviour of both local 
and nonlocal observables for integrable quenches in theo- 
ries with purely diagonal scattering, e.g. the sine-Gordon 
model at a reflectionless point, can be described by an 
appropriate GGE. 
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